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Abstract 

A latin bitrade (T°, T®) is a pair of partial latin squares which are 
disjoint, occupy the same set of non-empty cells, and whose correspond- 
ing rows and columns contain the same set of entries. A genus may be 
associated to a latin bitrade by constructing an embedding of the under- 
lying graph in an oriented surface. We report computational enumeration 
results on the number of spherical (genus 0) latin bitrades up to size 24. 

1 Introduction 

A latin bitrade (T°, T®) is a pair of partial latin squares which are disjoint, 
occupy the same set of non-empty cells, and whose corresponding rows and 
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columns contain the same set of entries. One of the earliest studies of latin 
bitrades appeared in [7], where they are referred to as exchangeable partial 
groupoids. Latin bitrades are prominent in the study of critical sets, which 
are minimal defining sets of latin squares ([2], [6], [16]) and the intersections be- 
tween latin squares ([11])- A genus may be associated to a latin bitrade, and 
those of genus zero are known as spherical latin bitrades. In this paper we report 
on the enumeration of the number of r-isomorphism classes of spherical latin 
bitrades up to size \T°\ = 24. 

For completeness we note that separated spherical latin bitrades are equiva- 
lent to cubic 3-connected bipartite planar graphs [5]. Our enumeration method 
does not consider equivalences where rows, columns, and symbols change their 
roles (this corresponds to changing the 3-colouring of the related cubic 3-connected 
bipartite graph), nor does it consider equivalences where (T°, T®) and (T®, T°) 
are equivalent (this corresponds to changing the 2-colouring of the related cubic 
3-connected bipartite graph). Thus our enumeration method counts a larger 
class of combinatorial objects compared to the method in [14]. 

2 Latin bitrades 

A partial latin square P of order n > is an n x n array where each e € 
{0, 1, . . . ,n — 1} appears at most once in each row, and at most once in each 
column. A latin square L of order n > is an n x n array where each e E 
{0, 1, . . . ,n — 1} appears exactly once in each row, and exactly once in each 
column. It is convenient to use setwise notation to refer to entries of a (partial) 
latin square, and we write k) E P if and only if symbol k appears in the 
intersection of row i and column j of P. In this manner, P C A\ x A 2 x A3 
for finite sets Ai, each of size n. It is also convenient to interpret a (partial) 
latin square as a multiplication table for a (partial) binary operator o, writing 
i o j = k if and only if (i, j, k) E T = T°. 

Definition 2.1. Let T°, T® C A x x A 2 x A 3 be two partial latin squares. Then 
(T°, T®) is a bitrade if the following three conditions are satisfied: 

(Ri) rnT® = 0; 

(R2) for all (ai, a 2 , 03) E T° and all r, s E {1, 2, 3}, r 7^ s, there exists a 
unique (61, b 2 , 63) G T® such that a r = b r and a s = b s ; 

(R3) for all (a 1; a 2 , a 3 ) E T® and all r, s E {1, 2, 3}, r ^ s, there exists a 
unique (61, b 2 , 63) e T° such that a r = 6 r and a s = b s . 

Conditions (R2) and (R3) imply that each row (column) of T° contains the 
same subset of A3 as the corresponding row (column) of T®. A bitrade (T°, T®) 
is indecomposable if whenever (J7°, J7®) is a bitrade such that U° C T° and 
[7® C T®, then (T°, T®) = (C/°, f/®). Bijections A4 -» for * = 1, 2, 3, give 
an isotopic bitrade, and permuting each Ai gives an autotopism. We refer to 
the bijections A x — > A 2 — > Ag, A3 — > A 3 as an isotopism. 



2 



In [9, 10] there is a representation of bitrades in terms of three permutations 
Tj acting on a finite set (see also [13] for another proof). For r G {1, 2, 3}, 
define the map f3 r : T® — > T° where (ai, a2, &3)/3 r = &2, ^3) if and only if 
a r 7^ fe r and a, = 6, for i ^ r. By Definition 2.1 each /3 r is a bijection. Then 
n, r 2 , r 3 : T° -> T° are defined by 

n = (3^/33, r 2 = /3 3 - 1 /3 1 , r 3 = /3r 1 /3 2 . (1) 

We refer to [n, T2, 73] as the Tj representation. We write Mov(7r) for the set of 
points that the (finite) permutation ir acts on. 

Generally we will assume that a bitrade is separated, that is, each row, 
column, and symbol is in bijection with a cycle of n, r 2 , and r 3 , respectively. 

Definition 2.2. Let t\, t 2 , T3 be (finite) permutations and let T = Mov(ti) U 
Mov(t 2 ) U Mov(t 3 ). Define four properties: 

(Tl) TlT 2 T 3 = 1; 

(T2) if pi is a cycle of n and is a cycle of tj then |Mov(/jj) n Mov(/9j)| < 1, 
for any 1 < i < j < 3; 

(T3) each r, is fixed-point-free; 

(T4) the group (n, T2, t 3 ) is transitive on T. 

By letting be the set of cycles of Tj, we have the following theorem, which 
relates Definition 2.1 and 2.2. 

Theorem 2.3 ([9]). A separated bitrade (T°, T* 8 ) is equivalent (up to iso- 
topism) to three permutations t\, t 2 , T3 acting on a set T satisfying (Tl), (T2), 
and (T3). // (T4) is also satisfied then the bitrade is indecomposable. 

To construct the Tj representation for a bitrade we simply evaluate Equa- 
tion (1). In the reverse direction we have the following construction: 

Construction 2.4 (r, to bitrade). Let t\, r 2 , r 3 be permutations satisfying 
Condition (Tl), (T2), and (T3). Let L = Mov(ti) U Mov(t 2 ) U Mov(t 3 ). Define 
Ai = {p I p is a cycle of t^} for i = 1, 2, 3. Now define two arrays T°, T* 81 : 

7 10 - {(pi, p 2 , As) A £ ^ and |Mov(pi) n Mov(p 2 ) n Mov(p 3 )| > 1} 
— {(pi, Pi-, Pi) I Pi € Ai and 2, a;', x" are distinct points of L such 
that xpi — x', x'p 2 — x" , x"ps = x}. 

By Theorem 2.3 (T°, T®) is a bitrade. 

Example 2.5. The smallest bitrade (T°, T* 8 ) is the intercalate, which has four 
entries. The bitrade is shown below: 
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The n representation is n = (000, 011)(101, 110), t 2 = (000, 101)(011, 110), 
r 3 = (000, 110) (Oil, 101), where we have written ijk for (i,j,k) G T° to make 
the presentation of the t* permutations clearer. By Construction 2.4 with T = 
{000,011, 101, 110} we can convert the Tj representation to a bitrade (J7°, t/®): 






(000,101) 


(011,110) 


= (000,011) 


(000,110) 


(011,101) 


(101,110) 


(011,101) 


(000,110) 


® 


(000,101) 


(011,110) 


= (000,011) 
(101,110) 


(011,101) 
(000,110) 


(000,110) 
(011,101) 



In this way we see that row of T° corresponds to row (000, 011) of t/°, which 
is the cycle (000, 011) of ti, and so on for the columns and symbols. 

Example 2.6. The following bitrade is spherical: 
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Here, the n representation is 

n = (000, 022, 044)(134, 142)(201, 213, 232, 220)(304, 333, 311) 
r 2 = (000,304,201)(213,311)(022,220)(134,232,333)(044,142) 
r 3 = (000, 220) (201, 311) (022, 232, 142) (213, 333) (044, 134, 304) 

If Z = [n, r 2 , t 3 ] then the inverse is denoted by Z~ 



T 2 TlJ 



This is equivalent to exchanging (T°, T ) for (T®, T°) and relabelling entries 
as the following lemma shows: 



Lemma 2.7. Let (T°, T®) 6e a separated bitrade with representation Z = 
[n, t 2 , t 3 ]. T/ien the inverse bitrade (T® , T°) /ias representation, denoted Z^ 1 . 
This representation is isomorphic to [r^ 1 , r^ 1 , r 2 Ti]. 

Proof. Observe that Definition 2.1 does not specify the order of the two partial 
latin squares in a bitrade so the bitrade (T®, T°) is well-defined. For this 
bitrade we have permutation representation [vi, ^ 2 , 1/3]. Using Equation (1) we 
find that V\ — fafi^ 1 and v 2 — PsPi 1 ■ Taking the conjugate of v r by fi^ 1 
gives us a representation that acts on the set T°. In particular, ft^vifls — 
/3s_ 1 p2p3 1 Ps =_ fe'fa =_ (^Ps)- 1 = and fc^lh = ^Ih^Pa = 
(3 1 1 /3 3 = (p3p-i 1 ) _1 = t 2 . This shows that the representation [1/1, i/ 2 , ^3] is 
isomorphic to [rf 1 , r^ 1 , t 2 ti]. □ 
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Figure 1: Rotation scheme for the combinatorial surface. Dashed lines represent 
one or more edges, where in is a cycle of n. The vertex in the centre is x. 

Let D be the graph with vertices T and directed edges (x, y) where XTi = y 
for some i. A rotation scheme can be imposed on this graph, by ordering 
the edges as shown in Figure 1. This turns D into an oriented combinatorial 
surface. Define order(ri) = z{t\) + z(j2) + z(t^), the total number of cycles, and 
size(ri) = |T|. By some basic counting arguments we find that there are size(ri) 
vertices, 3 • size(ri) edges, and order(r i ) + size(Tj) faces. Then Euler's formula 
V-E + F = 2-2g gives 



where g is the genus of the combinatorial surface. We say that the bitrade T 
has genus g. If g = then we say that the bitrade is spherical. 

3 Slide expansion 

Let [ti, T2, T3] be a spherical bitrade on the set T. Choose some ieT and fix 
a direction j e {1,2, 3}. Set k = j + 1 and I = k + 1 (calculating modulo 3). 
Suppose that (1) the j-cycle at x has length at least 3 and (2) the fc-cycle 
through x and the Z-cycle through w = xtj do not meet in any common point. 
Then the slide expansion at x € T in direction j is the following augmentation 
of the cycles of t\, t 2 , and r 3 by a new point u £ T: 



We denote the result of the slide expansion on X = [n , r 2 , T3] at vertex x in 
direction j by slide (X, x, j). 



orderfYi) = size(ri) + 2 — 2g 



(2) 



Tj : (a, x, w, b, ... ) (a, u, b, . . . ) and (x, w) 
T k : (...,x,z, ...) h-> (. ..,x,u,z, ...) 
Ti : (■■■,y,w,...) h-> (...,2/, u, «;,...)• 
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A separated bitrade [t±, T2, T3] is bicyclic if there exists j such that Tj con- 
sists of only two cycles. These bicyclic bitrades form the root nodes of our 
enumeration procedure due to the following result: 

Theorem 3.1 ([8]). Let Z = [n, T2, T3] be a spherical bitrade. Then there 
exists a sequence Wo, W\, Wg = Z where Wo is a bicyclic bitrade and 
Wi+i is the result of a slide expansion on either Wi or W[~ , for < i < I. 

The opposite of a slide expansion is a slide contraction. For a spherical 
bitrade [n, T2, t 3 ] on the set T, choose some u e T and fix a direction j e 
{1,2,3}. Set k = j + 1 and Z = k + 1 (calculating modulo 3). The slide 
contraction at u in direction j modifies the cycles of n, r 2 , and t 3 as follows: 

Tj : (a, u, b, . . . ) and (x, w) 1— » (a, x,w,b, . . .) 

T k : (...,x,u, z, ...) 

7"/: (■■■,y,u,w, ...) (...,2/, «;,...). 

A necessary condition for a slide contraction to be valid is that the w-cycle in 
Tfc and 77 has length at least 3 (this avoids the creation of a fixed point when 
u is removed). In general it is necessary to check the bitrade conditions (Tl), 
(T2), and (T3) to ensure that a slide contraction is valid. 

4 Enumeration 

To enumerate separated spherical bitrades we must fix the notion of r-isomorphism. 
Two bitrades [n, T2, T3] and [v\, f 2 , ^3] on the same set T are r-isomorphic if 
there is a permutation 9 € Sym(r) such that rf = Vi for all i e {1,2,3}. A 
r-automorphism of a bitrade is defined similarly. 

We follow the presentation in [15]. In the context of an algorithm for enu- 
merating combinatorial objects, the domain of a search is the finite set £1 that 
contains all objects considered by the search. The search space is conveniently 
modelled by a rooted tree, with nodes corresponding to elements of fi, and nodes 
joined by an edge if they arc related by one search step. The root node is the 
starting node of the search. We write C(X) for the set of all child nodes of X, 
and p(X) for the parent node of X. 

To reduce the search time we use the method of canonical augmentation [17], 
following the presentation in [15]. In terms of the search tree, the ordered pair 
(X,p(X)) characterises the augmentation used to generate X from the node 
p(X) during the search. In our case, this augmentation (X,p(X)) is the result 
of a slide expansion on p(X) . The goal of the canonical augmentation algorithm 
is to choose one of the possible parents of X to be the canonical parent . Formally, 
for any nonroot node X, the canonical parent m(X) e ft must satisfy: 

(CI) for all nonroot objects X, Y it holds that X = Y implies (A, m(X)) = 
(Y,m(Y)). 

The next property captures the canonical augmentations in the search tree: 



G 



(C2) a node Z occurring in the search tree is generated by canonical augmen- 
tation if (Z,p(Z)) ^ (Z,m(Z)). 

Algorithm 1 gives pseudo-code for the general canonical augmentation search 
algorithm. Algorithm 2 is a more explicit version for enumerating separated 
spherical latin bitrades. The correctness of Algorithm 1 relies on two properties: 

(11) for all nodes X, Y it holds that it X = Y, then for every Z G C(X) there 
exists aWe C(Y) such that (Z, X) ^ (W, Y). 

(12) for every nonroot node X, there exists a node Y such that X = Y and 
(Y,m(Y)) = (Y,p(Y)). 

Condition (II) says that isomorphic nodes have isomorphic children. Condi- 
tion (12) says that every nonroot node is generated by some canonical augmen- 
tation. 

Algorithm 1 CANAUG-TRAVERSE(A) 
1: report X (if applicable) 

2: for all Z G {C(X) n {aZ: a G Aut(A)}: Z G C(X)} do 

3: select any Z € Z 

4: if (Z,p(Z)) ^ (Z,m(Z)) then 

5: CANAUG-TRAVERSE(Z) 

6: end if 

7: end for 



Algorithm 2 CAN AUG- SPHERIC AL (A ) 
1: report X (if applicable) 

2: for all Z G {C(A) n {aZ: a G Aut(A)}: Z G C(X)} do 
3: select any Z € Z 
4: it (Z,p(Z)) = (Z,m(Z)) then 
5: CANAUG-SPHERICAL(Z) 
6: if has no parents then 
7: CANAUG-SPHERICAL^- 1 ) 
8: end if 
9: end if 
10: end for 

11: if (\X\ + 1) = (mod 2) then 

12: for all bicyclic Z with \Z\ = \X\ + 1 do 

13: CANAUG-SPHERICAL(Z) 

14: end for 

15: end if 



Theorem 4.1 ([15]). When implemented on a search tree satisfying (II) and 
(12), the algorithm CANAUG-TRAVERSE reports exactly one node from every 
isomorphism class of nodes. 
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4.1 Canonical form 

The canonical form Z is a relabelling of Z such that Z = Y if and only if 
Z = 1". We need to be able to efficiently compute the canonical form in order 
to construct the canonical parent of a bitrade. 

To compute the canonical form of a bitrade Z — [n, T2, T3], we fix a starting 
point x e T. We then perform a breadth-first traversal of the bitrade's under- 
lying directed graph. Vertices arc labelled from the set {1, . . . , \Z\} in the order 
that they are visited. At a vertex v we first visit the t\ cycle at v, then the 
t 2 cycle at v, and then the r 3 cycle at v (only if the cycle has not been visited 
before). As we traverse a cycle we label unseen vertices and append the cycle 
(using the new vertex labels) to an array C . After traversing a cycle we append 
a marker (here, —1) to denote the end of the cycle. While traversing a n cycle 
at the vertex v, we check if the neighbouring Tj and tu cycle at v has been 
visited (j, k =/= i). If it has not, we push (j, v) and (k, v) onto the queue, making 
sure that j < k. The lexicographically least C = C(x) is the canonical form Z 
of Z. Since Z encodes just the cycle structure of the n permutations, it follows 
that Z = Z' if and only if Z = Z' . 

The canonical parent m(Z) is now computed as follows. First we find the 
canonical form Z of Z. In our case, we choose the lexicographically maximal 
pair (j, u) such that a slide contraction may be performed on Z at u in direction 
j. Let Z be the result of this contraction. Then m(Z) = k~ 1 (Z ) where 
k: Z Z is the canonical relabelling map. 

4.2 Correctness 

Condition (II) holds since a slide expansion slide(X, x, j) is equivalent to slide(Y, y, j) 
where X = Y and y is the image of x under the isomorphism. In other words, 
the slide expansion relies only on the cycle structure of the bitrade, not its spe- 
cific labelling. Condition (12) holds because the canonical parent m(Z) is chosen 
only based on the canonical form of Z and so is independent of isomorphism. 

Algorithm 2 produces bicyclic bitrades at each step where even-sized bitrades 
are being generated. The following lemma and corollary are necessary to show 
that no r-isomorphism class of bicyclic bitrades is produced more than once. 

Lemma 4.2 (Lemma 10.1 [8]). Let [t\, r 2 , r 3 ] be a bitrade such that Tj, for 
some j G {1,2,3}, consists of exactly two cycles. Let one of the cycles be 
(xo, . . . ,x n -i). Then the other cycle of Tj can be expressed as (x' n _ 1 , . . . , x' ) 
in such a way that Tj+i consists of cycles (2^,2:^) and Tj-\ consists of cycles 
(Xi +1 ,x't), for0<i<n. 

Corollary 4.3. Let Z be a bicyclic bitrade. Then there is no X such that 
slide(A, x, j) = Z or slide(A, x, j) = Z^ 1 for any x, j. 

Proof. The slide contraction shortens the cycle at some point u in r^+i and 
Tj_i. By Lemma 4.2 at least one of the cycles in Tj +1 or Tj_i has length 2, so 
the slide contraction would introduce a fixed point, contradicting (T3). □ 
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It can be shown that for any bicyclic bitrade Z we have Z = Z 1 . Corol- 
lary 4.3 shows that lines 11-15 of Algorithm 2 are correct. 

5 Implementation and relation to other work 

See [12] for C++ computer code. The MPI library is used to allow parallel 
processing of the search tree on a Linux cluster. The size of the automorphism 
group Aut(Z) is bounded by the size of T. Further, the number of homogeneous 
bitrades is quite small, so usually the automorphism group is even smaller. We 
just compute the elements of Aut(Z) explicitly. 

Cavenagh and Lisonek [5] showed that spherical bitrades are equivalent to 
planar Eulerian triangulations. To verify our results (for small n) we used 
plantri [3] to generate planar Eulerian triangulations. We then applied a si- 
multaneous vertex 3-colouring and face 2-colouring, giving a bitrade (T°, T®). 
Then each conjugate of (T°, T®) was produced, and finally isomorphic copies 
removed. 

The number of r-isomorphism classes of spherical bitrades with \T°\ = n is 
given below: 
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The third column is the value e^{n v ), which is the number of isomorphism 
classes of eulerian plane triangulations with connectivity at least 3 (see Table 6 
in [4] for values n/ = 2n in their first column and also [1]). The enumera- 
tion resulting in e^(n v ) considers isomorphic any two graphs that are related 
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by a permutation of the face 2-colouring, changing the orientation of the per- 
mutations, and permuting the cyclic ordering. So for large values of n there 
will be at most 2 • 2 • 3 = 12 ways to relabel one of their graphs to obtain 
one of our bitrades. For example, with n = 24 we have ez{n v ) — 24873248 
and 12 • 24873248 = 298478976 w 297753519. Also, Wanless [18] enumerated 
spherical bitrades up to size 19, under various various equivalences, but not 
T-isomorphism as studied in this paper. 
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